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This  paper  presents  a  procedure  for  the  design  of 
superconducting  electrical  machines.   The  magneto- 
static  field  problem  of  a  cylindrical  superconducting 
machine  is  solved.   Equations  for  the  magnetic  fields 
are  developed,  and  from  these,  expressions  for  the 
various  machine  parameters  are  obtained.   These  expres- 
sions are  adapted  to  a  computer  solution  of  the  design 
problem.   The  computer  program  obtains  the  design 
parameters  for  a  minimum  volume  machine  design. 

A  design  study  for  a  typical  marine  electrical 
propulsion  system  is  conducted.   The  results  of  this 
study  indicate  quantitatively  the  reduction  in  weight 
end  space  to  be  obtained  in  using  superconducting 
electrical  machines  vice  conventional  electrical 
machinery.   The  study  al3o  presents  some  general 
characteristics  of  superconducting  machines. 
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Figure  1    Machine  cross  section 
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I.   INTRODUCTION 

A.  Background 

The  relative  flexibility  and  simplicity  accruing  to 
electrical  power  systems,  in  comparison  to  mechanical  power 
systems,  has  long  been  recognized.   Specifically,  electrical 
power,  and  transmission  systems  provide  flexibility  of 
installation,  are  easily  suited  to  automation  and  have  a  high 
degree  of  dependability. 

Despite  these  advantages,  and  others  that  could  be  obtained, 
the  full  benefit  of  electrical  power  systems  has  not  been  realized 
in  marine  propulsion  applications.  To  be  sure,  electric  drive 
systems  have  been  used  in  ship  propulsion,  however,  such  use  has 
been  mostly  limited  to  small  or  intermediate  size  propulsion 
plants.  The  primary  reason  that  electric  drive  systems  have  not 
been  adopted  for  large  poller  plants  has  been  that  the  weight  and 
space  for  such  systems  is  greater  than  that  required  for  geared 
propulsion  systems.   In  addition,  increased  cost  (as  much  as  20$) 
and  lower  efficiency  (85$  for  DC)  have  been  contributory  factors 

It  is  advantageous  to  consider  some  of  the  benefits  that 
could  be  obtained  with  electric  propulsion  in  order  that  aspects 
of  its  utilization  might  be  judged  in  proper  perspective. 

As  mentioned  above,  an  electric  propulsion  system  affords 
greater  flexibility  of  installation.   In  constrast  with  a 
conventional  geared  propulsion  system,  extensive  shafting  is  not 

'Numbers  in  parantheses  indicate  applicable  reference  in  the 
bibliography. 
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required  and  the  location  of  the  power  source  is  not  as 
restricted.   In  addition,  the  large,  noisy  reduction  gear  could 
be  eliminated.   The  noise  generated  by  this  component  is  of 

increasing  military  significance. 

.v. 

If  a  cycloconvertor  is  usod,  additional  benefits  derive. 

It  would  be  possible  to  use  a  constant  frequency  generator, 
thus  permitting  a  more  efficient  design  of  the  prime-mover. 
This  will  also  permit  greater  use  of  newly  developed  power 
sources  such  as  the  gas  turbine.   Added  to  this  would  be  more 

adaptation  of  remote  and  automatic  controls,  plus  much  improved 

(2  9) 

maneuvering  characteristics.   ,7/ 

Thus  it  seems  that  if  the  disadvantages  (weight,  space,  and 
possibly  efficiency)  of  electric  propulsion  are  reduced,  it 
will  be  given  greater  consideration  as  an  alternative  to  geared 
propulsion  in  larger  power  systems. 

^ •  The  Role  of  Superconductors 

It  has  been  realized  for  some  time  that  when  the  temperature 
of  an  electrical  conductor  is  reduced  tc  a  few  degrees  absolute, 
the  resistance  of  the  conductor  decroascs  to  an  unmeasurable 
amount.  The  implications  of  this  fact  have  considerable  value 
for  use  in  electrical  machines.   If  the  resistance  in  the  windings 
of  an  electrical  generator  or  motor  is  removed,  much  greater 
current  densities  can  be  achieved  with  an  accompanying  increase 
in  the  magnetic  flux  densities  within  the  machine.  ->l      If  this 


"An  electronic  power  trans  .  .-ion  device  that  affords  chang   in 
voltage  and  frequency.   See  reference  [{.. 
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incroase  in  flux  density  is  obtained,  then  greater  power  can 
be  generated  within  a  given  volume.  Thus  it  follows  that  the 

volume  of  an  electric  machine  can  be  reduced  if  superconducting 

(3) 

windings  are  used.   Studies  to  date,  ^ '   have  demonstrated  the 

significant  reduction  in  volume  and  weight  in  superconducting 
machines. 

However,  the  application  of  superconductivity  is  not  ifithout 
its  complications.   The  first  question  which  must  be  asked 
concerns  the  method  of  operating  a  machine  wherein  all,  or  at 
least  some,  of  its  parts  are  at  a  temperature  in  the  region  of 
liquid  helium. 

It  is  at  this  point  convenient  to  limit  the  scope  of  the 
examination.  The  ensuing  discussion  will  treat  only  machines 
with  a  rotating  field  winding,  operating  in  the  superconducting 
region,  enclosed  by  a  stationary  armature  at  ambient  temperature. 
This  particular  configuration  is  only  one  of  several  that  could 
be  considered. 

There  are  two  particular  reasons  why  this  configuration  is 
selected.   The  more  obvious  is  that  it  would  be  impractical  to 
transfer  the  high  machine  power  output  through  slip  rings,  as  would 
be  necessary  if  the  armature  were  rotating.  The  other  reason 
results  from  the  fact  that  large  heat  losses  are  encountered 
when  superconductors  are  subjected  bo  alternating  currents.  *•*' 
Thus,  only  the  field  winding  is  made  superconducting. 

Nonetheless,  most  of  the  techniques  in  this  paper  are 
generally  applicable  to  any  cylindric*       ;rical  device. 
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With  the  discussion  constrained  to  considering  a  rotating 
superconducting  field  winding,  it  is  possible  to  dilineate  some 
design  procedure.   Thus  far,  it  has  been  found  that  the  use  of 
a  thermal  bottle,  or  dewar" ,  filled  with  liquid  helium  and 
encompassing  the  rotor,  is  a  satisfactory  solution  to  the 
problem  of  maintaining  the  windings  in  a  superconducting  state. 
Several  techniques  are  under  development  for  introducing  the 
helium  into  the  rotating  dewar  and  passing  electrical  leads 
into  the  low- temperature  region.   Generally  it  is  favorable 
to  use  the  liquid  helium  dewar  in  conjunction  with  a  liquid 
nitrogen  dei/ar  immediately  exterior.  The  helium  maintains  the 
rotor  at  about  I}..20K  while  the  nitrogen  offers  a  77°K  buffer 
for  improved  cryogenic  performance.   At  first,  one  might  think 
that  the  refrigeration  load  for  maintaining  such  low  temperatures 
would  be  excessive.  However,  it  should  be  realized  that  once 
the  low  temperatures  are  obtained,  there  is  very  little  heat 
generated  by  the  windings.   Indeed,  heat  introduced  into  the 
low-temperature  region  throtigh  the  shaft  and  support  structure 
is  more  significant.   A  reasonable  estimate  of  the  power 
requirements  for  refrigeration  during  operation  of  a  super- 
conducting  machine  would  bo  about  0.5  percent  of  rated  power. 

In  addition  to  the  requirement  that  a  superconductor  be 
maintained  at  a  low  temperature,  there  :1s  a  limitation  on  the 
amount  of  current  carried  while  exposed  to  a  magnetic  field. 
As  indicated  in  Figure  2,  if  the  magnetic  flux  density  increases, 


'Named  after  the  inventor. 
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the  permissible  current  density  is  reduced.   If  current  or 
field  exceed  the  limits  imposed  (the  so-called  critical 
values),  the  conductor  reverts  to  its  normal  state.   Thus 
it  is  necessary,  in  any  particular  machine  design,  to 
determine  that  combination  of  maximum  flux  density  and 
current  density  which  best  suit  the  design. 


flu*  density 


Zil2i£iL-i-   Current -Flux  Relationship  for  a 
Typical  Superconductor 


The  very  fact  that  a  superconducting  machine  develops 
extremely  intense  magnetic  fields  invites  yet  another 
complication.   Since  fields  of  the  order*  of  several  kilo- 
gauss  can  be  expected  within  the  machine,  disturbance  of 
equipment  in  the  vicinity  of  the  machine  will  certainly 
result,  as  well  as  unbalanced  loads  on  the  superconducting 
field  winding.   Therefore,  it  is  necessary  that  a  shield  be 
provided  to  confine  the  magnetic  fields  within  the  machine. 
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This  shield  would  typically  be  of  laminated  steel  con- 
struction and  will  be  a  significant  portion  of  the  size 
and  weight  of  the  machine. 

V •      Design  Program 

Though  superconducting  machines  of  the  size  suggested 
by  the  foregoing  discussion  have  not  yet  been  built,  they 
are  certainly  feasible.   Anticipating  the  construction  of 
such  machines,  it  is  considered  desirable  to  investigate 
some  of  their  expected  properties.   Since  there  is  vested 
interest  in  the  reduced  size  and  ;</eight  of  these  machines, 
it  is  natural  that  a  design  procedure  be  developed  for 
obtaining  minimum  weight  and  volume. 

The  objective  herein  is  to  find  a  procedure  for  finding 
design  characteristics  of  least  size  machines.   The  procedure 
developed  is  founded  on  a  field-theory  solution  of  the 
electric  machine  and  is  adapted  to  a  computer  solution. 
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II .   PROCEDURE 

A.   Development  of  a  Model 

The  first  step  in  the  analysis  was  to  obtain  a 
mathematical  model  for  the  magnetic  fields  involved.   This 
model  was  necessary  to  determine  the  parameters  which  would 
be  used  in  machine  design.   As  mentioned  previously,  the  model 
and  attendant  equations  are  for  a  rotating  field  winding  inside 
a  three-phase  stator  winding. 

Analysis  proceeded  by  considering  both  the  field  and  stator 
windings  to  consist  of  cylindrical  current  sheets.  Equations 
for  the  magnetic  fields  generated  by  each  winding  were  developed 
separately  and  later  combined  to  obtain  necessary  inductance  and 
energy  expressions. 

To  illustrate  the  analysis,  consider  the  field  windings'. 
As  illustrated  in  Figure  3*  a  cylindrical  current  sheet,  K(r). 
of  thickness  dp  ,  located  at  some  radius  r,  R.jZrZR0,  is  examined 


Figure  3  Distribution  of  Current  Sheet  K(r) 
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Fourier  analysis  of  this  cur-rent  sheet  yields: 


£ 


L^dr 
K(r)  =  [__,       -~—     sinnp0,  (II-l) 

n  odd   nTr 


Knowing  that  these  current  sheets  exist  beWeen  r  a  R.  and 
r  =  Rp,  and  that  the  fields  associated  with  each  sheet  arc 
Laplacian,  enables  determination  of  the  field  expressions.  >y 
After  finding  an  expression  for  the  magnetic  field  for  r  -  Rp, 
one  can  find  the  contribution  due  to  the  shield.   Since  at 
r  =  R  ,  with  an  infinitely  permeable  shield, 


ns0      "f0  U±-^ 


and  at  r  -  0, 


H    -  finite  (II-3) 


an  expression  is  found  for  tine  fields  induced  by  the  shield. 
The  expression  thus  obtained  is  then  added  to  that  for  the 
fields  generated  by  the  field  winding. 

This  same  procedure  obtains  for  determining  the  magnel  '• 
field  system  resulting  from  the  armature  windings ,  The  only 
difference  involved  is  in  the  current  sheet  model  used,  since 
the  armature  is  s  three-phase  winding,  with  60-degree  phs 
bolts. 
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3p 


Figure  l\.         Distribution  of  Armature  Current  Sheet 


Fourier  analysis  of  the  armature  current  sheet  yields 


Ka(r) 


n  odd 


kJ  dr 

•a  nrr        .  ^ 

— — —     cos  -«-     sin  np0   . 


(II-14-) 


Expressions  for  magnetic  field  systems  of  the  field  and 
armature  windings  are  presented  in  Tables  1  and  2. 

Using  the  equations  for  the  magnetic  fields t    an  expression 
for  energy  stored  in  the  armature  field  system  is  obtained  from 


(2x 


W 


a 


fRs 


I  F-o   (H*   +  Ilf0  )  lp  dp  d0  .     (II-5) 
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Equating  the  result  in 


wnere 


Wft  =|Lalf,  (11-6) 


I  -  ff(Ro  ~  *1}  Ja  ,  (II-7) 

a       g_     _ 

a 


yields  an  expression  for  the  armature  self -inductance.   A 
similar  calculation  is  performed  to  obtain  the  field  winding 
self- inductance . 

Mutual  inductance  is  found  from 

IfLaf  "^  '  (II~8) 


where  A  is  the  total  field-winding  flux  linkage  of  the  armature 

windings.   By  recognizing  that  6N  /rr(R   -  Rt)  expresses  the 

a    o     j. 

number  of  armature  turns  per  unit  cross-section,  then  the 
second  spatial  differential  of  the  flux  linkage  is  expressed 

by 

?  6N 

d^A   =   — %  a        Mop  dp   d0  H~     lap   df.  (11-9) 

rr(Pv     -  Rf )  x  y  tt        -P 

0  1  0 
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TABL?J  2 


Magnetic  Flo Id  Eauations  for 


Phase  A  of  Armature  Windings 
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i    R^   -   R? 
1     -!•    -r 


*  R|   ln(R0/Ri) 


p  cos   20 


p  sin  20 


Hrt   = 


iirt   = 


Ld 


Tf 


4 


1    - 


Hi' 


Ro  Rn   ~   Ri 

-*•   In  -2     +  .  < 


R, 


■  o4 


4R 


4R' 


pcos    20 


psin  20 


R0<P<Kr 


Hp  - 


l0 


i-rx  ^uo 


4) 


^s 


ivS 


cos    20 


sin  20 


-15- 

Integrating  this  expression  over  the  cross-section  of  the 
armature  yields  the  total  flux  linkage  and  thus  the  mutual 
inductance.   Expressions  for  field  and  armature  3elf- 
inductance  and  mutual  inductance  are  presented  in  Table  3. 

B.   Design  Equations 

If  one  considers  the  amplitudes  of  the  line-to-neutral 
armature  voltage  and  armature  current,  one  obtains 

Va  =WeMIf,  (11-10) 


T    rr(R2  -  R?)J 

*„  =     °    i   a 

a   — gF- 


(11-11) 


where 


P  =  |v  I  , 
g   2  a  a' 


1       2Nj 


(II-12) 


(II-13) 


This  power  is  internally  generated  power,  which  must 
be  reduced,  to  account  for  armature  impedance  and  the  load 
power  factor,  in  order  to  obtain  rated  output  power. 


Figure  5>   Vector  Diagram  of  Generator  Voltages 
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U  =  -= 


16>U0NJ1 


773P(P2-*)d-y2)2 


La  - 


16^0N|1 


TT3P(P2-^)(1-X2) 


TABLE  3 

Inductance  Equations 
P  t   2 

(p-2)    -    (p+2)y4+   4yP+2 
+   2^(1-yP+2) 

(p-2)         (p*2)x4+  4xp+2 


'Hp|2p 


-prl(-^2)(g) 


7r3p(l~x2)(1~y2)         Ih0J 


1-x-P+2   +    1-xP*2    / RA  2p 


P   =   2 


Lf   - 


LQ    = 


8  /40N|l 


irO-y^r 


TT30-x2)2 


4  >    I  3-2  I  2P 


l(l.yA)   +y4lny  +  1(1.r)^ 


i(1-/)+Anx+i(,-xV(| 


2p 


6/i0NaNfl(Uy2)    fug1 
TT3(1-x2)  [Ro 


ml  +  iizadi  (52 

x  A  lag. 


2r, 
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From  the  vector  diagram  in  Figure   5>>   where 


xs  =    <Vv 


(Il-llf.) 


the  ratio  V. /V  is  obtained.   Then  rated  power  is 
"C  a 


p  =  p  (v. /v  )  cos  e. 

g  u  a 


(II-15) 


In  developing  the  equations  describing  the  machine,  some 
approximation  was  necessary  to  account  for  the  effects  of  the 
end- turns.    '   Figure  6  shows  the  geometry  involved  in  making 
this  approximation. 


Fi£ure_^6   End-Turn  Geometry 


If  one  considers  an  outer  armature  turn  at  the  mean  radius  of 
a  given  phase  winding,  it  is  soon  th    he  angular  displacemen 
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of  this   coil   is  l|ir/3p  radians.      Since   this   coil   is   at   a  radius 


r  -  |(RC+    R±)  (11-16) 


the  approximation  can  be  made  that 


Using  this  expression,  it  has  been  found*   that  a  reasonable 
length  to  assume  for  self- inductance  calculations  is  the  total 
length 

lt  =  1  +  2A1,  (11-18) 


whereas  the  length  to  be  used  for  mutual  inductance  is 


1  =  1+  Al.  (11-19) 

m 


As  described  in  the  Introduction,  the  low  resistance 
character  of  a  superconducting  coil  is  not  only  a  function  of 
temperature,  but  also  depends  on  its  current  density  and  the 
magnetic  field  to  which  it  is  exposed.   Thus  for  a  given  design 
it  is  necessary  to  determine  the  maxi   .•  magnetic  flux  density 
associated  with  any  particular  current  density,  to  ensure  that 
the  superconducting  field  winding  will  not  revert  to  its  normal 
state.   The  particular  mechanics  employed  to  observe  this 
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limitation  are  detailed  in  the  program  construction  techniques 

of  Appendix  C . 

As  also  mentioned  in  the  Introduction,  the  machine  design 

considered  employs  a  shield  (constructed  of  laminated  steel) 

which  contains  the  dense  magnetic  fields  within  the  machine. 

For  this  shield  to  be  effective,  the  fields  it  constrains  must 

be  less  than  those  which  will  cause  saturation  of  the  shield. 

If  the  shield  material  saturates  at  some  density  B   .  , ,  then 

*      rated 

one  limitation  is 


B  „  {r  -  R  )  -  b   .  ,. 
max. v     o     rated 


(11-20) 


Another  limitation  exists  in  the  requirement  that  the  maximum 
flux  density  carried  by  the  shield  must  never  exceed  B   .  ,, 
As  illustrated  in 


Figure  7  Shield  Flux  Pattern 
Figure  7»  the  total  flux  accumulated  by  the  shield  over  an  angle 
of  ff/2p  must  not  be  sufficient  to  cause  saturation  in  the  cross- 
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section  at  A-A.      Since 


frr/2p 
®A-A=|C  VV1Rod*'  (II"21) 


a  final  constraint  is  obtained  in 


®A-A  -  Brated  <HoS  "  W  <"-22> 


There  are  two  final  considerations  which  may  govern  a 
machine  design.   First,  the  heat  dissipation  in  the  armature 
winding  must  be  considered.   Generally  the  maximum  permissible 
current  will  be  determined  from  the  limitations  imposed  by  the 
method  adopted  for  cooling  the  armaturo.  For  purposes  of  this 
design  it  was  assumed  that  an  average  armature  current  density 
of  1.6  x  10  amp/m  could  be  accepted  with  chilled  water  cooling 
of  the  tending.   Finally,  one  must  consider  the  speed  of  the 
machine  rotating  parts.  The  strength  of  the  r.otor  structure 
will  provide  a  limit  to  the  rotor  tip  speed  that  can  be  sustained, 
On  consideration  of  previous  conventional  machine  design  it  was 
determined  that  700  ft/sec  was  a  reasonable  upper  value  for  tip 
speed.  However,  since  methods  of  constructing  superconducting 
windings  are  far  from  being  fixed,  this  value  for  maximum  tip 
speed  is  speculative.   The  tip  speed  of  the  rotor  will  also  bo 
governed  by  the  performance  of  the  rotating  dewar.   In  this 
respect,  the  fluidic  and  thermodynamic  properties  of  liquid 
helium,  under  high  centrifugal  forces,  will  govern.      Since  the 
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aspects  of  this  consideration  are  yet  to  be  determined,  it  is 
assumed  here  that  the  previously  stated  maximum  tip  speed  is 
acceptable. 
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III.   RESULTS 

The  significance  of  this  work  lies  in  the  results  of  a 
design  study  which  was  conducted,  employing  the  design 
procedure  described  in  Part  II.   The  design  study  involved 
what  was  considered  a  typical  marine  electrical  propulsion 
system;  a  turbine  driven  generator  rated  at  26  MW  (35*000  hp) 
at  3600  rpm,  and  two  direct  driva  propulsion  motors  rated  at 


17 > 500  blip  (13  MW)  at  300  rpm.   Other  parameters  assumed  were: 


J .«    =      1.0  x  10   amp/m  ' 


J    =      1.6  x  106  amp/m2 
a 

Bmax  =      5.0  web/rc2 

BratecT     2"°  w9b/»2 
p^         -        lj.ooo  kg/irr 

p  =  l\.000  kg/ra^ 

(copper  with   ,\\$  packing  factor) 

ps  8000  kg/in3 

(steel) 

S  ~  .  05>m 

Tip  speeds  of  up  to  21I4.  m/sec  (700  ft/sec)  were  considered 
in  bo tli  the  generator  and  motor  designs.   Also,  for  the 
generator,  designs  with  po\;er  factors  of  1.0..  .8^0  and  .707 
were  examined.   Design  characteristics  of  machines  with 
1  to  6  pole  pair:  were  computed  for  all  variations. 
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The  computer  results  and  graphs  of  the  characteristics 
of  some  of  the  more  interesting  designs  are  included  in 
the  following  pages. 

One  quantity,  not  calculated  by  the  computer  program, 
that  should  be  examined  is  the  machine  efficiency. 

Since 

R  «  P/aH  (III-l) 

a   -g— 

where  P  =  resistivity"  1.9  x  10"*'  ohm-m  and  cross-section 
area  is 

A  =  «d-*2)Ro  *  -!±   ,  (III-2) 

x   

a 
for  a  .1|5  packing  factor,  then  for  the  optimum  generator 
design  (pf  =  .707,  P  =  6,  V..   =  211+  m/sec) 

R  /H2  -  6.7  x  10"8  ohm/turn2, 
a  a 

Knowing  K  I  ,  it  is  possible  to  find 
a  a 


A   =  965  watts. 


a  a 
These  losses,  however,  are  for  only  one  phase  and  were 

computed  with  a  peak  value  for  I   instead  of  rcns .   Therefore, 

a 

resistive  losses  for  all  three  phases  are 

I2R  :  965  '  |  r':  lli50  watts. 
a  a         C- 

If  an  equal  amount  is  allowed  for  eddy  current  losses,  a 

figure  of  2900  watts  is  obtained  for  the  total  losses. 

Since  this  is  for  a  26  KW  generator,  this  loss  will  cause 

a  negligible  reduction  in  efficiency. 
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Figuro  9 
Motor   Specific  Weight 
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IV.   DISCUSSION  OF  RESULTS 

The  most  significant  result  of  this  study  lies  in  the 
low  values  of  specific  weight  (lb/hp)  for  superconducting 
machines.   As  indicated  in  Figure  8  specific  weights  on  the 
order  of  0.25  lb/hp"  are  theoretically  possible  in  machines 
of  the  size  typically  found  in  marine  propulsion  plants. 

The  results  of  the  study  of  direct-drive  propulsion 
motors  were  less  conclusive.   It  is  apparent  from  Figure  9 
that  specific  weights  of  I.I4.O  lb/hp  are  possible  and  that 
much  lower  specific  weights  are  likely  to  be  achieved  in 
machines  with  larger  numbers  of  poles  and  higher  tip  speeds. 
Nonetheless  a  specific  weight  of  l.Ij.O  is  still  an  improvement 
over  a  conventional  machine.  This  is  evident  from  a 
comparison  with  the  data  for  the  propulsion  motor  in  the 
USS  Hunlcy,  AS-31,  given  in  Table  l\.. 

hp  15,000 
pole  pairs  26 
rpm  If?  7 

lb/hp  10.9 

Table  ij   Characteristics  of  Propulsion  Kotor 
of  USS  Ilunley  AS-31*11* 

Consideration  of  the  results  of  the  generator  designs 
indicates  that  in  high  speed  machines,  maximum  tip  spoed  will 


"Compared  to  l\. . 0  lb/hp  for-  conventiona]  generators. 
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be  a  limit  in  machine  design.  However,  this  is  not  likely 

to  be  the  case  in  lower  speed  machines  such  as  the  direct 

drive  motor. 

Machine  characteristics  such  as  low  x  and  high  efficiency 

a 

are  expected  as  a  result  of  the  smaller  amount  of  armature 
winding  required. 

Aside  from  the  significance  of  the  data  given  for  the 
various  machine  designs,  there  were  tx«ro  other  notable 
results  of  this  work.   As  descr5.bed  in  the  procedure  for 
determining  the  field  winding  design  (detailed  in  Appendix  B), 
in  order  to  reduce  the  number  of  variables  involved  it  wa3 
assumed  that  the  fields  induced  by  the  shield  did  not  affect 
B     .   However,  the  study  conducted  in  Appendix  E  indicates 
that  in  some  instances  (especially  for  small  values  of  p) 
this  assumption  is  not  warranted.   From  Figure  16  it  is 
found  that  for  the  optimum  generator  design  (x  =  .822, 
y  =  .8)4.8,  p  =  6)  the  error  introduced  is  about  one  percent, 
whereas  for  the  optimum  motor  (x  =  .801,  y  =  .'887,  p  =  6) 
the  error  is  about  .80  percent.   These  accuracies  are  ex- 
pected to  be  much  greater  than  the  accuracy  of  the  math- 
ematical model  being  used. 

Another  result  that  became  evident  during  the  design 

study  was  that  the  linear  interpolation  technique  used  to 

find  the  minimum  volume  design  (described  in  Appendix  D) 

x>ras  not  sufficiently  accurate  in  all  instances.   Since  the 

minima  of  the  volume  vs ,  R  curve  is  quite  broad  for 

o 
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p  =  1  end  2,  the  linear  interpolation  yielded  erroneous 
results  in  a  few  designs. 

The  accuracy  of  solution  provided  by  the  design  program 
(within  the  limits  imposed  by  the  assumptions  as  discussed 
in  Appendices  B  and  E)  is  arbitrary.   In  the  design  study, 
a  variation  of  1%   in  output  power  was  allowed.   Likewise, 
in  instances  of  solving  for  maximum  flux  density,  the  rate 
of  change  was  required  to  be  less  than  10"  *"  at  the  calculated 
maximum.   By  changing  these  tolerences,  the  designer  can 
obtain  the  accuracy  desired  in  the  solution. 


-35- 

V.   CONCLUSIONS  AND  RECOKKENDATIONS 

From  the  results  of  this  study,  the  following  can  be 
concluded: 

1.  Superconducting  electrical  machines  should  provide 
considerable  savings  of  weight  and  space  in  electrical 
propulsion  systems. 

2.  The  electrical  efficiency  of  superconducting  machines 
is  substantially  higher  than  the  efficiency  of  conventional 
machines. 

3.  The  solution  of  the  design  problem  for  a  super- 
conducting machine  can  bo  made  at  least  as  accurate  as  the 
design  of  a  conventional  machine, 

ij..  High  speed  machines  will  be  limited  in  speed  of 
rotation  by  the  maximum  tip  speed  that  can  be  sustained  by 
the  rotating  parts. 

5.   An  optimum  design  for  a  slow  speed  machine  will 
require  a  large  number  of  poles. 

In  the  conduct  of  further  investigation  of  the  optimal 
design  problem,  or  in  application  of  the  design  program 
presented  here,  the  following  is  recommended: 

1.  The  body  of  data  (array  B(p,  y) )  representing  the 
graph  of  B   .   /const,  vs.  y,  should  be  expanded  to  account 
for  larger  values  of  p. 

2.  The  design  program  should  be  modified  to  allow 
variation  in  VTIP.   As  indicated  by  Figures  8  and  9j  an 
optimum  machine  design  may  occur  for  some  value  of  VTJP 
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that  is  less  than  the  maximum  VTIP.   By  searching  over 

a  range  of  this  variable,  a  numerical  technique  could 

easily  be  used  to  find  the  optimum  design. 

3.   When  a  design  is  obtained,  it  is  necessary  to 

ensure  that  the  error  in  B_    ,  as  described  in 

p  max 

Figures  13  to  16,  is  not  excessive.   It  may  be  desirable 
to  incorporate  these  figures  in  the  design  program,  as 
was  done  with  array  B(p,  y).   It  would  then  be  possible 
to  check  each  design  within  the  program,  making  appropriate 
adjustment  in  the  design  parameters  until  a  satisfactory 
result  is  obtained. 

1|.  To  ensure  greater  accuracy  in  the  subroutine  Mill, 

(7) 
described  in  Appendix  D,  a  false  position  solution  technique 

should  be  substituted  for  the  linear  interpolation  method. 

A  basic  routine,  such  as  the  one  presented  on  the  following 

page,  might  be  considered. 
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APPENDIX  A 

Input  Data 

In  organizing  the  machine  program  that  was  ultimately 
used,  it  was  necessary  to  anticipate  what  information  was 
likely  to  be  available  or  fixed  for  the  machine  design.   It 
was  anticipated  that  certainly  the  desired  power  output 
would  be  known.   The  angular  velocity  of  the  rotor  and 
power  factor  of  the  load  may  also  be  determined  from 
particulars  of  the  intended  machine  application.   Further, 
it  was  expected  that  materials  selection  x;ould  precede  this 
portion  of  a  design,  therefore  various  material  character- 
istics will  be  available.   Also  expected  is  some  knowledge 
of  construction  techniques  to  be  employed,  e.g.  the  type  of 
armature  cooling  used  and  characteristics  of  the  rotor 
design.   Finally,  some  information  concerning  the  dcwa.r  is 
necessary.   Specifically  the  anticipated  thickness  of  the 
dewar  wall  must  be  provided.   Thus,  input  data  will  include: 

a.  power  output,  P 

b.  mechanical  frequency,  COm 

c.  maximum  armature  current  density,  Jr 

a 

d.  dowar  wall  thickness,  A 

e.  maximum  tip  speed  of  dewar,  V. . 

(fixed  by  strength  considerations  in  rotor  design  or 
limitations  imposed  by  dewar  technology) 

f.  power  factor,  pf 


-fcO- 

g.   flux  density  to  cause  shield  saturation,  Brate(^ 

h.   average  density  of  rotor,  armature  and  shield,  p  „,  p  ,  n 

As  mentioned  in  the  Introduction,  the  superconducting 
material  used  in  the  field  windings  will  have  characteristic 
values  of  critical  current  and  critical  flux  density.  There- 
fore, additional  input  data  will  be  a  value  of  current  density 
in  the  field  windings,  Jf,  and  the  associated  value  of  critical 

flux  density,  B 

max 

For  any  field  winding  design,  the  geometry  of  the  magnetic 

fields  can  be  found.   If  J^,  is  known,  the  maximum  flux  density 

can  bo  calculated.  This  maximum  flux  density  must  be  found 

to  ensure  that  it  does  not  exceed  the  critical  flux  density 

associated  with  J~   (B     ) . 

■*•    m  a  x 

Conversely,  if  J^,  and  B  m  8  x  are  given,  certain 
limitations  on  field  winding  design  are  imposed. 

It  was  found  expedient  to  provide  a  body  of  data  which 

could  be  used  in  the  computer  program  for  finding  the  best 

field  winding  design  for  some  given  J„  and  B     „   The 

i       ma  x 

method  for  obtaining  and  utilizing  this  data  is  described 
in  Appendix  B. 


-in- 


APPENDIX  B 


Input  Data  for  Field  Winding  Design 

In  the  design  of  the  field  windings,  three  assumptions 
are  made.   First  it  is  assumed  that  the  fundamental  component 
of  the  magnetic  fields  is  the  significant  part.   This  assump- 
tion is  fairly  standard  and  has  been  shown  to  be  acceptable. 
Second,  it  is  assumed  that  the  fields  generated  by  the  arma- 
ture may  be  neglected.   Since  the  armature  fields  oppose  the 
fields  generated  by  the  field  windings,  this  approximation 
is  conservative.   Third,  it  is  assumed  that  the  fields 
reflected  by  the  shield  may  be  neglected.   Although  these 
fields  add  to  those  from  the  field  winding,  they  are  small, 
becoming  insignificant  as  p  increases.   These  approximations 
are  made  in  order  to  decrease  the  number  of  variables  in  the 
calculations . 

From  Table  I,  it  is  seen  that  in  the  region  R-.  LqL   R? 
the  maximum  field  intensity  will  be  in  the  p-rdirection 
for   0-0.   In  accordance  with  the  above  approximations, 


u 


for  p  r  2, 


He  = 


or 


B 

const. 


2Jf 
rr(^-p2) 


BpTT 

2^0JfR2 

r>  ■  1 1 


-2P  •■  (2  •  p) 


(^-p:2-(2-P)P>-2 


R. 

P/ 
(B-l) 


Ui~P^) 


2p  |   +(2-p) 
2 


J 


R 


2 


P-l  p-i-2 

-(2-p)y   X 


(B~2) 


.A 


-lj.2- 


If ,  for  a  given  value  of  y, 


d(Bp  /const. ) 

dTy7Rp 


0 


(B-3) 


is   solved  for    (  P/R~),    and  the  value   obtained  used  to 
evaluate  Eq.    (B-2),    the  maximum  value   of  Bp  /const,    is 
obtained  for   that  y. 


For  the   case  p   -   2 


and 


H      = 


B 


23. 


TT 


R 


k 

b1^ 


BpTf 


const . 


VA 


R. 


In 


Q 


If    R0     7 


+  p 

R. 


In 


(B-l^) 


(B-5) 


Thus,    if  R„   is  known,    it   is   possible   to   obtain 


B 


P  max , 


f(p»y). 


(B-6) 


In  the  following  computer  program,  for  the  IBM  1130; 

B     /const,  is  commuted  for  values  of  p  =  2  to  6  and 
p  max 

y  =  0.0£  to  0. 95 »  The  results  are  displayed  graphically 

in  Figure  10.   Using  these  results,  it  is  possible  to  find 

the  value  of  y  corresponding  to  a  given  3  _     /const., 
J  ta      &       p  max      ! ' 

thus  fixing  the  radial  dimensions  of  the  field  windings. 
The  case  p  -  1  is  not  computed  since  for  this  value  of  p 


E 


/const.  -  ].. 


(B-7) 


-to- 

In  the  program,  the  following  notation  is  used: 

p  ■  number  of  pole  pairs 

y  =  R -,/Rpj  R-,  is  the  inside  field  radius,  etc. 

R  =  P  /Rp,  i.e.  Rl  in  the  program  is  the  first  value 
of  p/Rp  and  is  not  the  inside  field  radius. 

B  =  Bp  /const, 

B  MAX  =  B  ^     /const. 
p  max 

The  program  procedure  is  simply  a  solution  of  Eq.  (B~3) 

and  evaluation  of  Eq,  (B-2)  or  Eq.  (B-5)  for  the  maximum 

(7) 

value.   For  the  case  p  ~  2  a  false  position  method  '   is 

used  to  solve  Eq.  (3-3) >  whereas  the  Kewton-Raphson 

/ON 

technique^  '  is  employed  for  all  cases  p  f.   2. 

The  printed  results  of  this  program  include  values  of 
p,  y  and  B^    /const,  which  are  used  in  the  design  program, 
The  data  represents  the  curves  in  Figure  10.   Within  the 
design  program,  B     /const,  is  computed,  with  which  it  is 
then  possible  to  obtain  from  the  curves  a  valve  of  y  which 
will  satisfy  the  constraint  imposed  by  the  current  and 
critical  flux  density. 

The  program  used  to  solve  for  the  data  in  the  array 
B(p>y)»  ~nd  the  subroutine  used  to  interpret  B(p,y)  to 
find  a  value  for  y,  are  included  in  the  following  pages. 
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APPENDIX  C 

Design  Program 

The  design  computer  program,  written  for  the  IBM  1130, 
is  an  iterative  solution  of  the  basic  design  equations  to 
find  that  machine  design  which  requires  the  least  volume. 

The  program  uses  the  normalized  variables  x  and  y  which 

represent  the  thickness  of  the  armature  and  field  windings. 

The  variables  R„  and  R   serve  to  fix  the  radial  location 

d  o 

of  the  windings.   The  machine  length  is  found  from  consider- 
ation of  the  working  volume  required  to  obtain  the  necessary 
internally  generated  power.   The  reduction  from  internal 
power  to  rated  power  is  obtained  from  Eq.  (11-15).  The 
iterative  procedure  followed  is  based  on  varying  values 

of  R  .  which  was  found  to  be  the  most  convenient  parameter 
O'  ^ 

for  this  purpose. 
Since 

the  variation  of  R   with  R  was  first  examined.   By 

OS         o  J 

evaluating  Eq.  (11-21)  and  substituting  in  Eq.  (11-22)  one 


(C-2) 


obtains 
R 

= 

W*(1-»p+2) 

fRi 

^oJ 

P 

1-2 
R2 

0 

+R 

ot> 

pP                         o 
*o 

;p) 

O 

C-3) 


-50- 


Figure  11   Variation  of  R   with  R 


For  snail  values  of  R  ,  the  shield  is  in  the  vicinity 

o 

of  intense  magnetic  fields  and  R   must  be  largo.   A3  R 
°  os  °        o 

increases,  the  thickness  of  the  shield  decreases  but  R 

1  os 

is  still  larger  than  R  .  The  value  of  R  which  requires 

the  smallest  R   is  found  to  be 
os 


R    ,.   =  (pC.  ) 
o  opt.    y*   1' 


I/d+1 


(c-k) 


Though  this  value  of  R  does  not  necessarily  yield  a 
minimum  volume  machine  (since  no  variation  of  1.  is 
considered),  it  is  generally  near  the  minimum  volume  value 
because  of  the  significance  of  R   in  Eq.  (C-l). 

03 

An  R  obtained  from  Eq.  (C-k)  will  satisfy-  the  constraint 
o 

of  Eq.  (11-22).  However,  it  must  also  satisfy  Eq.  (11-20). 


-Sl- 


Since 


B  v(r=R  ) 


.  w. 


T.'(2*p) 


/R  1  p+2 

,-    p+2i    "2 

\  0» 


o  —  rated 


(c-5) 


it  is  required  that 


Ro  i    (C1/p)1/P+1i 


(C-6) 


Since 


O  Opt 


=  PCn  *  -1 


(c-7) 


Eq.  (11-20)  is  also  satisfied. 

Turning  now  to  the  design  program,  the  following 
procedure  is  carried  out  for  successive  valties  of  p: 

From  the  input  data,  one  obtains 

R.   -  V.  .  /a)    , 
i      tip  wm' 


(C~8) 


R, 


=  R  -  A 

1    *- 


(c-9) 


Evaluating 


B  ,  ,  it 

max      „    2  a  J  ..R0 

"o  i  2 


(C--10) 


and  entering  the  curves  of  B     /const,  vs.  y,  described 
^  max  J 

in  Appendix  B,  a  value  of  y  is  obtained,   tfith  this 

information  it  is  possible  to  evaluate  C.  and  find  R    . 

^  1  o  opt . 

which  is  used  to  calculate  an  initial  machine  design, 


Having  selected  a  value  for  R   one  obtains 

o 

R. 
-      l 

X  -  R~  • 


(C-ll) 


At  this  point  in  the  program,  the  difference  between  the 
internally  generated  power  and  the  desired  power  output  is 
unknown.   Assuming  that  this  difference  is  largely  due  to 
load  power  factor,  it  is  approximated  that 


P„  =  P  /pf 
g     o 


(C-12) 


Since 


xvhere 


5   =  i   V  I 
g    2  a  a 

S  I  «eHVa 


(0-13) 

(C-lij.) 


We  =  Pwm 

i*  =ff(1-y2)Rl  Jf 

2N,, 


(c-i5) 
(0-16) 


=*(l-x2)R2ja  (0-17) 

a 
and  the  appropriate  expression  for  M  is  selected  from 
Table  II,  one  obtains 

f  f?  \P+2 


P   =6^„Vf   (l-yP+2)  f^l 


^  x 


1-x 
li~p 


-P42 


(2+p)2 


K 


(C-18) 


for  p  J*  2,  or 


-53- 


L 


In  ± 

x 


1   +  l-x~- 


T 


(C-19) 

for  p  =  2.   From  Eq.  (C-l8)  or  Eq.  (C-19),  it  is  possible  to 

solve  for  1  .  Then,  in  accordance  with  Ecu  (II-18)  and 
m 

Eq.  (11-19), 


1.  =  1+  Al 

t    m 


=  1 


TT 


"m 


3p  <1+X>V 


(C-20) 
(C-21) 


It  is  now  possible  to  find  the  reduction  in  power  from 
P  to  P  .   Applying  the  law  of  cosines  to  the  diagram  in 
Figure  5>  results  in 


Normalizing,  one  obtains 


2  ,  ..2 


(v./v  r  \   x  -  2  v./v  x  cos  (5  f  e)  -  i  =  o 

t   a      a      t   a  a      2 


(C-23) 


where 


x  =  x  I  /V 
a    s   a  a 


(C-2W 


is  the  per  unit  armature  reactance  loss.   If  Eq.  (C-23) 

is  solved  for  V./V  ,  then  Ec  .  (II-1E>)  can  be  evaluated 

t,      a 

for  P  . 
o 

The  value  of  P  thus  obtained  is  compared  with  the 
o 

value  specified  in  the  input  data.   If  the  difference  is 
greater  than  a  given  tolerance  (1$  of   the  desired  value, 


for  bhis  program),  the  value  of  P  assumed  in  Jtnq «  (C- 

S 


-A- 

is  increased  (or  decreased)  by  the  percentage  difference 
between  actual  and  desired  output  power  (see  statement  30 
in  the  design  program).   The  program  then  returns  to  Sq. 
(C-l8)  or  Sq.  (C-19),  and  calculation  continues  in  this 
fashion  until  an  acceptable  difference  is  reached. 

VJhen  a  satisfactory  P  is  obtained  the  program  continues 
with  calculation  of 


°1 

R  a   =  --  +  R 
os   Rp      ( 


(c-3) 


Vol.  =  TT  R_   1. 
OS     t 


(C-l) 


Weight  =  Trlt((l-y2)R|  pf  +  (l-x2)n2  9&   +(R2s-R2)  Pg ) 


V,/N   =  V^/V      '    V  /N   =  V./V 


a 


a 


(G-25) 


am  //pi^d^  2)Rf2  Jf 

TT2(1-X2)R] 


X 


i>.«p 


for  p  i  2,    or 


V./N   =  V./V 
t   a    -c   a 


(in  I 

X 


1-x 
(2+p); 


p+2 


TT2(l-x2)R2 


f-  1  ,„   Llm 

5  u-x  )} 


(C-26) 


(0-27) 


-SS- 

for  p  =  2,  and 

tt(1-x2)R2  J 

HI  *  °  a  .  (C-28) 

a  a        6 

The  calculations  described  above  are  carried  out  for  the 
initial  value  of  R  .  The  program  continues  by  incrementing 
R   (by  0.1  in  this  program),  returning  to  Eq.  (C-ll)  and 
repeating  all  the  steps  through  Eq.  (C-28).   R  continues 
to  be  incremented  until  it  is  found  that  the  volume  for 
successive  designs  is  increasing,  which  indicates  that  the 
minimum  volume  design  is  within  the  range  of  those  designs 
that  have  been  completed.   The  program  then  calls  a  sub- 
routine, described  in  Appendix  D,  which  evaluates  the 
data  for  the  designs  that  have  been  corslet ed,  and  from 
this  data  determines  the  minimum  volume  design. 

Within  the  design  program,  the  following  notation  is  used 

PWR  =  desired  P 

o 

POUT  =  actual  P 

o 

B(p>y)  -  array  of  data  representing  the  graphs 
in  Figure  10 

OUTPT  b  array  in  which  parameters  of  ench  design 
are  stored 

BC  =  B      ,    . 
£>  max/const. 

XA  =  un-normalized  armature  reactance,  which 
becomes  normalised  when  stored  in  OUTPT. 

VA01T  =  V  /N 
a   a 

VTON  =  V.  /H 
t   a 


-56- 


PG  = 

:  P 
g 

ROS 

s  R 

OS 

VOL 

"  Volume 

WGT 

=  Uoight 

HIA 

=  N  I 
a  a 

A  print -out  of  the  source  deck  for  the  design  program 
is  included  in  the  following  pages. 
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APPENDIX  D 

Subroutine  to  Find  the  Minimum  Volurae  Design 

After*  the  design  program  computes  three  designs  of 
successively  increasing  volume,  it  calls  the  subroutine 
MIN.   The  data  stored,  in  OUTPT  at  this  time  will  include 
at  least  two  designs  beyond  the  minimum  volume  design. 


Vol . 


^Opt   ^01   ^01        Ro3 


04 


Figure  12a   Data  Stored  in  OUTPT 


dVOL 
dR0 


VUL1  DVQL2-— 


DV0L3 


Calculated  by  I  '. 
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The  program  computes  dVol/dR   for  the  three  values  of 

R  that  bracket  the  minimum  volume  design, 
o 

Since 

Vo1  =  *  Ros    h  (D-1J 

then 

dVol   .  „  R2   dlt    +  w.   d(Ros)  ,n  ^ 

o  o  o 

Considering  Equations  (C-l8),  (C-19),  and  (C-21),  it  is 

possible  to  express 

p 

1   =  &     rP"2  +  *    (r   +h.)  (D-3) 

t   C     o     3p    o    i'  \      si 

inhere 

Co  =  f(x).  (D-lj.) 

Therefore 

£  =  (D„2)  ~S  rP-3  _   g o  Rp-2  +  tt_     ,  e) 

o  o        C     o  ■/- 

o 

and ,  from  Eq .  ( C - 3 ) > 

:   -2pC2R~2p'"1  h  2C-(l-t))R""p  4-  2R  .    (D-6) 

J.  O  ±    "    0        o 


d(R2J 

OS' 


dR 
o 

When  values  of  DVOL  for  R^,  ,  Rrt0  and  R„0  have  been 

computed,  linear  interpolation  is  used  to  find  the  value 

of  R  for  which  DVOL  (R  )  =0.   Control  then  returns  to 
o  o 

the  main  design  program  x^iere  the  design  characteristics 

for  th:i  s  v.     of  R   ore  found, 

o 
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A  print-out  of  tho  source  deck  for  this  sub-routine  is 
included  in  the  follx^oing  pages. 
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APPENDIX  E 

Effect  of  the  Shield  on  Maximum  ..Flux  Density 

As  described  in  Appendix  3,  the  assumption  was  made, 
in  determining  characteristics  of  the  field  windings,  that 
the  effect  of  the  shielding  could  be  neglected.   The 
significance  of  this  assumption  is  examined  as  follows: 

If  B p  /const,  is  computed,  as  in  Appendix  B,  without 
neglecting  shield  effect,  the  result  is 
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Eq,  (E-l)  end.   Eq.  (E-2)  can  bo  written  in  terms  cf  the 
aalized  parameters  x,  y  and  o  = 
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Assuming  that  0  is  constant  over  a  range  of  machine 

sizes,  it  is  possible  to  compute  the  actual  B  ^    /const. 

'       r  p  max 

in  terms  of  x  and  y.  The  solution  technique  is  the  same 

as  that  described  in  Appendix  B.   One  can  then  compare 

the  value  of  B      /const,  obtained  from  the  simplified 

p  max 

equations  with  that  value  obtained  from  the  method 
described  above. 

The  results  of  such  a  comparison  are  displayed  in  the 
following  graphs.  The  data  is  for  the  cases  of  p  =  3  to  6, 
which  are  of  most  interest. 
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Figure    13 
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Figure  15 
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Figure  16 
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